
A comparison between high-energy radiation 
background models and SPENVIS trapped- 
particle radiation models 





GOAL: 

Determine the background 
rate due to cosmogenic and 
terrestrial background 
radiation for the next 
generation of X-ray and 
Cosmic-Ray missions. 

SPENVIS use is highlighted. 
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The Problem 





• The next generation of space-based X-ray and Cosmic-ray experiments are designed 
for unprecedented signal sensitivity — > high background rejection. 

• The magnitude of the background must be determined as well as mitigation 
strategies developed, both in the choice of instrument design and signal processing 
capability. 

• Simulation packages (GEANT) can be used to accurately determine the impact of 
background radiation if the spectra of all backgrounds are known (and the Fermi 
team has developed an extensive set of background models (Mizuno et al., ApJ 614)). 

• This problem is more acute for the new X-ray experiments, historically it has been 
difficult to simulate the radiation background that mimic X-ray signals with GEANT and 
the low-energy background in low-Earth orbit, e.g. ISS, is more uncertain. 

• For mission planning, SPENVIS has proved and invaluable tool to determine the orbit 
fraction at specific geomagnetic latitude bands and determine where the trapped- 
particle flux is too large to allow physics operation — > helps define duty cycle. 

• For X-ray missions, better knowledge of the low-energy (electron and proton) 
terrestrial background is needed ... can SPENVIS help? 
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Background Flux Modeling 


Instrument Modeling 
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• Experiments in orbit are irradiated by multiple sources: 

• Cosmic Ray Primaries: protons, helium, e 1 , heavy nuclei. 

• Diffuse X-ray and Gamma-ray background. 


• Atmospheric 'Backsplash' and trapped particles: protons, photons, e 1 , neutrons. 

• High values of trapped-particle flux preclude science operation. 

• Fluxes and spectra depend on orbit parameters (altitude and geomagnetic latitude) and 
solar activity (Solar Modulation Factor) 

• Fluxes and spectra sampled by radiation-sensitive detectors also depends on surrounding 
material, which acts as both a shield and a target. 

• Particle physics based Monte Carlo simulations provide mechanism to assess impact of 
radiation background on an experimental measurement. 
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The Neutron star Interior Composition ExploreR 

(NICER) 





Front-side hotspot rotates through the line of sight 


Simplified MGEANT model 


NICER is an array of 56 individual telescopes, 
each with grazing-incident optics and a silicon 
drift detector with a ~1.1 m separation. 

X-ray signal energy range: 200 eV - 10 keV 

NICER will measure radii of neutron stars, 
which reveals the interior composition. 
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Increasing compactness (M/R) and light bending 



High-Fidelity Mass Models 




Developed by Steve Sturner @ GSFC 

• MGEANT mass model based 
upon CAD instrument model. 

• Materials in model match that 
in instrument, i.e. alloys vs pure 
Al (X-ray fluorescence). 

• Care taken to model detector 
housing & materials. 


Window (Si 3 N 4 + Al) Mu | ti _ Layer internal Collimator (Al, Ti, Cr, W) 

Inactive Detector Material (Si) ^ Capsule (Ni) 


tector Substrate (Alumina) 
Getter (ST172 + Mo) 



Base Plate (Au plated Kovar) 


Pins (Au plated Kovar) 


Stud (Kovar) 
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li, (cm 1 ) 


MGEANT Modeling 

We are using MGEANT, which is based upon GEANT3, but with improvements: 

• MGEANT uses GLECS: an extension program for the GEANT3 simulation package that 
incorporates improved physics into the simulation of Compton and Rayleigh scattering. 

• In GLECS, the total cross sections for an element are interpolated from the Evaluated 
Photon Data Library (EPDL [2]) tables. These data include electron binding effects and are 
valid from 10 eV to 100 GeV. 




• MGEANT used forTGRS on WIND, GRIS, and SPI on Integral. 

• "We hove successfully used MGEANT with the lower energy bounds set as low as 200 eV." 
-MGEANT Manual. 


• For GEMS 'A atm DME polarimeter, comparison to GEANT4 simulations and laboratory 
data show excellent agreement (NIM paper currently being written). 



Energy (key) 


Comparison of silicon linear attenuation coefficient (dots) to NIST values. 
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Curves: NIST XCOM values 
Red Points: a-sampling > 200 eV 
Blue Points: a-sampling > 1 keV 

Ml = P <*Tot /( uA )/ U ~ 1-66 X 10 ' 24 g 



Background Radiation Models 



• based upon Mizuno et al., ApJ 614 (2004) & Fermi updates and are data driven. 

• includes dependence on geomagnetic latitude, extended to lower energy for some. 

• considers background sources outside of regions of high fluxes of trapped particles. 





Energy (GeV) Energy (MeV) 


• these spectra can then be used in GEANT to generate background event files. 

• the definition of the event generating surface and the integral flux defines an event rate. 

• event selection criteria can then be applied to minimizing the background rate while 
maximizing the signal acceptance, requires a instrument response model. 
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d“(|)/dEdQ (particles/(m MeV s sr)) 



Secondary Proton Model 





Mizuno Downward Secondary Proton Flux (0 < 0 M < 0.9) 
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2ndory Proton Envelope 

»/ <t> = 1 .31 E-2 protons/ (cm 2 s sr) 


Mizuno201 1 (0 < 0 W < 0.2) 

/ 4 = 7.34E-3 protons/(cm 2 s sr) 

Mizuno201 1 (0.2 < 0 M < 0.3) 
/ 4> = 3.53E-3 protons/(cm 2 s sr) 

Mizuno201 1 (0.3 < 0 M < 0.4) 
/ 4 = 2.95E-3 protons/(cm 2 s sr) 

Mizuno201 1 (0.4 < 0 M < 0.5) 
/ 4 = 2.67E-3 protons/ (cm 2 s sr) 


Mizuno201 1 (0.6 < 0 M < 0.7) 

/ 4 = 3.46E-3 protons/ (cm 2 s sr) 

Mizuno201 1 (0.7 < 0 M < 0.8) 

/ 4 = 4.94E-3 protons/ (cm 2 s sr) 

Mizuno201 1 (0.8 < 0 M < 0.9) 

/ 4 = 8.80E-3 protons/(cm 2 s sr) 

Mizuno CR Protons (®J= 
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Kinetic Energy (GeV) 


Secondary Proton Flux is a function of 
geomagnetic location (& up * down at high 6 M ) 

• Integral flux goes from 2.7 — > 8.8 (rel units). 

• Spectral shapes rather similar. 

• Define envelope function to throw MC events, 
then need to scale by relative flux in appropriate 
bands of geomagnetic latitude. 

MIZUNO ET AL. 

TABLE 1 


Model Parameters for Proton Secondary Flux at Satellite Altitude, Based on AMS Data 


Region 

Direction 

Model 

Parameters 3 

o.o<e M < 0.2 

Downward/ Upward 

Cutoff PL 

0.136/0.123/0.155/0.51 

0.2 < 0 M < 0.3 

Downward/ Upward 

Broken PL 

0.1/0.87/600/2.53 

0.3 <6 U < 0.4 

Downward/ Upward 

Broken PL 

0.1/1.09/600/2.40 

0.4 <0 M < 0.5 

Downward/ Upward 

Broken PL 

0.1/1.19/600/2.54 

0.5 <0 M < 0.6 

Downward/ Upward 

Broken PL 

0.1/1.18/400/2.31 

0.6 <0 M < 0.7 

Downward 

Broken PL 

0.13/1.1/300/2.25 


Upward 


0.13/1.1/300/2.95 

0.7 < 6m< 0.8 

Downward 

Broken PL 

0.2/1.5/400/1.85 


Upward 


0.2/1.5/400/4.16 

0.8 <0 M < 0.9 

Downward 

Cutoff PL 

0.23/0.017/1.83/0.177 


Upward 

Broken PL 

0.23/1.53/400/4.68 

0.9 <0m< 1.0 

Downward 

Cutoff PL 

0.44/0.037/1.98/0.21 


Upward 

Broken PL 

0.44/2.25/400/3.09 


a Fo/a/£bk(GeV)/ft for the broken PL (PL) model (eq. [7]) and Fo /F\ /a/£ c (GeV ) for the cutoff PL 
model (eq. [8]). F 0 and F\ are given in counts s _l m~ 2 sr -1 MeV 1 . 
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SPENVIS use I: 

Defining geomagnetic orbit fraction 





• Using the '< Coordinate generators/Spacecraft 
trajectories' tool the orbit position parameters 
(latitude & longitude) can be generated in equal 
time spacing based on altitude (400 km), 
inclination (51.6°), and duration (30 days). 

• Assuming a location of the Earth's magnetic north 
pole (X p =85.9,(|)p =147.0) the geomagnetic latitude 
(A) for a given orbit (X,4>) can be calculated: 

sin A = sinAsinA F +cosAcosA p cos (0 - 0 F ) 


The geomagnetic latitude fractions can be 
thus calculated and be used to determine the 
background rate vs geomagnetic latitude. 


GeoMagnetic 
Latitude Bin 
(radians) 

Fraction of 
Orbit within 
Bin 

Running Total 

0 - 0.1 

0.082 

0.082 

0.1 -0.2 

0.083 

0.165 

0.2 -0.3 

0.084 

0.249 

0.3 -0.4 

0.086 

0.335 

0.4 -0.5 

0.089 

0.424 

0.5 -0.6 

0.095 

0.519 

0.6 -0.7 

0.107 

0.626 

0.7 -0.8 

0.134 

0.76 

0.8 -0.9 

0.173 

0.933 

0.9 -1.0 

0.067 

1 


ISS Orbit fraction versus 
Geomagnetic Latitude obtained 
from Spenvis Analysis 
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SPENVIS use II: 





Excluding regions of high flux of trapped particles 


• Using the ' Radiation sources and 
effects /Trapped proton and electron 
fluxes' tool, the integral flux vs energy 
of trapped electron (AE8) and protons 
(AP8) can get generated as a function 
of the orbit position parameters 
(latitude & longitude). 

• Selection of a maximum integral flux 
vs position can then be used to 
exclude regions of the orbit where 
trapped flux is too high to allow 
science operation — > duty cycle. 

• Note: ROSAT (53° inc, 580 km) was 
able to operate in upper high flux 
region. 



Longitude 


AE8 (SolarMax) map of integral electron 
flux (E Elec > 40 keV) for 1 day ISS orbit. 
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Latitude (Degrees) 



Duty Cycle Determination 





Using the results of the NICER background modeling (outside of regions of 
high trapped particle: electrons > 1 Hz/cm 2 above 40 keV; protons > 1 Hz/cm 2 
above 100 keV) and excluding the trapped particle flux regions — > duty cycle. 



Legend 

Green: backgound < 
0.01 cts/s/keV 

66% of ISS orbit 

Blue: background > 
0.01 cts/s/keV 

17% of ISS orbit 

Red: trapped particle 
flux > 1 Hz/cm 2 

17% of ISS orbit 
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Background Rate vs GeoMag Lat 





GeoMagnetic 
Latitude Bin 
(radians) 

Fraction of 
Orbit within 
Bin 

Orbit Running 
Total 

Background Rate within 
Bin (CPS) 

(0.2 keV < E dep < 2 keV) 
Goal: < 0.02 CPS 

Background Rate within 
Bin (CPS) 

(2 keV < E dep < 10 keV) 
Goal: < 0.08 CPS 

O 

1 

O 

0.082 

0.082 

1.35E-2 

4.98E-2 

O 

i-> 

1 

o 

k> 

0.083 

0.165 

1.36E-2 

5.01E-2 

0.2 -0.3 

0.084 

0.249 

1.37E-2 

5.03E-2 

o 

1 

cn 

o 

0.086 

0.335 

1.39E-2 

5.14E-2 

0 

1 

o 

In 

0.089 

0.424 

1.45E-2 

5.41E-2 

0.5 -0.6 

0.095 

0.519 

1.45E-2 

5.40E-2 

0.6 -0.7 

0.107 

0.626 

1.70E-2 

6.48E-2 

0.7 -0.8 

0.134 

0.760 

1.83E-2 

7.15E-2 

0.8 -0.9 

0.173 

0.933 

2.14E-2 

8.59E-2 

O 

T— 1 
1 

CD 

o 

0.067 

1 

2.76E-2 

1.14E-1 
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d“(|)/dEdQ (particles/(m MeV s sr)) 


SPENVIS Use III: Low-energy flux guidance 




Fermi-based background flux models extrapolation 
to low-energies (Fermi interested in > 10 MeV). 


Radiation Source 

Energy Range 

Cosmic Ray Primaries (p, a, e~, e + ) 

100 MeV- 200 GeV 

Cosmogenic Photons 

1 keV — 100 GeV 

Atmospheric Secondaries (p, n, e~, e + ) 

10 MeV -10 GeV 

Atmospheric Photons 

30 keV-IOOGeV 


Need to extend flux models to lower energies, 
where input spectrum x instrument response 

— > 0 signal (analysis shows this is true for high 
energies and is mainly true at lower energies, 
but would like to extend further down in E). 

Electron CSDA range in Al: 

100 keV: 70 um; 1 MeV: 2 mm; 10 MeV: 2.2 cm 
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Integral Flux (cm 2 a -1 ) 


SPENVIS Flux Reports 



AE-S MAX Orbit averaged flux 



0.01 0.10 1.00 10.00 

Energy (MeV) 


ISS Orbit averaged differential and 
integral electron flux. Leads to 
over-prediction of background rate. 


• SPENVIS models go to lower E: 

• 40 keV for electrons. 

• 100 keV for protons. 

•To assess low-energy contribution to 
background, need to get differential 
flux as a function of orbit position, 
where trapped radiation flux is falling 
to 0. 

• SPENVIS reports integral electron 
and proton trapped flux as a function 
of position, but LSB is 1 Hz/cm 2 
assumed isotropic 

— > 8e-2 particles/(s sr cm 2 ) 

• Outside areas of trapped particles 
integral flux range (for ISS): 

< IE-2 - 2E-1 particles/(s sr cm 2 ) 
assuming E MIN ~ 10 MeV. 
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Flux Extrapolation to Low-E 



Atmospheric 'secondary' e + e' 

• Outside regions of high flux of 
trapped particles. 

• Spectra defined by AMS-01 
measurements with E > 100 MeV, 
but with extrapolation to 10 MeV. 

• Extrapolation to 10 keV yields wide 
range of flux values. 

• Background rate prediction thus 
vary over a wide range. 


10 10 

May 22, 2013 


★Estimate scaling low flux bit of 
AE8max model of 1 Hz/cm 2 in a 
0.04 - 0.1 MeV energy bin. 

— > Naive extrapolations to lower 
E overpredicts the background. 
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Summary 


• We have been using GEANT with high-fidelity mass models and robust background flux 
models to predict the background rate in the next generation of X-ray & Cosmic-ray 
experiments: the flux models and procedures are robust. 

• For X-ray mission modeling, we would like to extend our background models to lower 
energy (down to ~100 keV), and SPENVIS offers the opportunity to do this extrapolation. 

• Future work: 

• Understand the fidelity of the SPENVIS trapped particle radiation models: 

• Are the SPENVIS trapped flux models of sufficient accuracy to provide flux predictions at 
low-flux edge of the trapped radiation? 

• As presented at this workshop, the AE9 and AP9 updates provide differential spectra 
and reports the trapped radiation flux to a factor of 10 lower than the AE8 and AP8 
models. This offers the opportunity to to use the position dependent flux predictions for 
the trapped particle models set low energy flux limits at the low trapped particle flux 
boundaries defined by the lowest non-zero reported flux for the trapped radiation. 
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CALET: Calorimetric Electron Telescope 
(New Cosmic Ray ISS Mission) 





Main Telescope: Calorimeter (CAL) 

• Electrons: 1 GeV - 20 TeV 

• Gamma-rays: 10 GeV - 10* TeV 
(Gamma-ray Bursts: > 1 GeV) 

• Protons and Heavy Ions: 

10’s of GeV - 1 ,000* TeV 

• Ultra Heavy (Z>28) nuclei: 

E > 600 MeV/nucleon 
Gamma-ray Burst Monitor (CGBM) 

• X-rays/Soft Gamma-rays: 7keV - 
20MeV 




CHD: 2 x,y layers of 45 cm long scintillator paddles. 
IMC: 16 x,y layers of scintillating fibers interleaved 
with tungsten plates (3 X 0 ) 

TASC: 12 x,y layers of lead tungstate logs (27 X 0 ) 
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